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Abstract 

Atoms trapped in optical lattice have long been a system of interest in the AMO community, and 
in recent years much study has been devoted to both short- and long-range coherence in this system, 
as well as to its possible applications to quantum information processing. Here we demonstrate for 
the first time complete determination of the quantum phase space distributions for an ensemble of 
^^Rh atoms in such a lattice, including a negative Wigner function for atoms in an inverted state. 
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I. INTRODUCTION 



Much of the excitement and promise of new fields such as quantum information processing 
would not be possible without the development of sophisticated techniques for manipulating 
and measuring quantum systems. In systems such as ion traps, quantum dots, Bose-Einstein 
condensates, and entangled photons, it is the fact that quantum states can be accurately 
prepared and observed which has enabled a wide range of experimental advances, and held 
out hope for the creation of practical quantum technologies. In such systems, quantum states 
inevitably evolve into partially mixed states, which can be characterized either by a density 
matrix, in the case of a finite-dimensional Hilbert space, or by phase space quasi-probability 
distributions such as the Wigner function[l|| or the Husimi distribution^. Techniques for 
extracting these functions, generally referred to as "quantum tomography," have long been^a 
topic of active research, in the context of the electromagnetic field! 
„e..a, .o^m, dissocating .oleCesfl. e.a„,led photonsQ 



Rydberg states!^ 
llJ, lUj, and ion trapsjjj 



One system which has led to many interesting effects and proposals is the "optical lattice," 
in which atoms are trapped in a periodic potential formed by a standing wave of light 
jjl^. Iinl. 11^ . This analog condensed-matter system has been used to study Bragg 



beams 



scattering[l2|, precision measurement using atom interferometryj3|, the fragmentation of 
Bose condensates 0|i squeezed states of atomic motion 2111. quantum feedback|22l|. the 
Mott insulator transition j2j|, and quantum logic gates|24. 3|, to name only a few 
examples. Typical probes of the system include the Bragg scattering of probe beams, which 
is sensitive to the localization (X^) of the atoms; the time-varying transmission of the 
lattice beams themselves, which is sensitive to the instantaneous force exerted on the atoms, 
related to the centre-of-mass position (X) in each well; and atom-interferometric probes of 
long-range coherence. Here we present for the first time a complete characterisation of the 
quantum state of atoms trapped in optical-lattice wells, by extracting the Husimi and Wigner 
distributions for atoms in several different initial states. In particular, we demonstrate the 
"nonclassical" Wigner function for atoms with a population inversion, whose negative value 
at the origin is analogous to that for the 1-photon Fock state and for the excited state of 
the single ion in a trap. 
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II. THEORY 



The Husimi distribution is a well known quasi-probability distribution related to a state's 
density matrix p by the equation Q [a] = ^ {a\p\a), where |a) is a coherent state, that 
is, a Gaussian distribution centered at a, with minimum-uncertainty width which can be 
defined by a specified harmonic oscillator. It has the practical advantage that its value for a 
particular a is an observable and can be obtained directly with one measurement. Because it 
is always real and always positive it is sometimes referred to as a classical quasi-probability 
distribution. 

The Wigner distribution, which has a one-to-one correspondence with the Husimi dis- 
tribution, is particularly useful for identifying nonclassical states such as Fock or inverted 
states, for which it takes on negative value (as in fact it does for any non-Gaussian states). It 
can be defined by the expression W {x,p) = ^ {x + q\ p\x — q) exp [—2ipq] dq 

0- The 

Wigner distribution is unique in that it allows one to determine the marginal probability 
distribution of either coordinate by integrating over the other: P {x) = J W {x, p) dp and 

P{p) = {x,p) dx. 

The measurement of these quasi-probability distributions becomes particularly simple in 
a harmonic oscillator. The Husimi distribution can be measured by evaluating the expression 

Q{\a\,e) = -{a\p\a) = -{0\Dma\)R {9) pR^ {6) D {\a\) |0) , (1) 

71 TT 

where 9 = arg {a}. That is, instead of projecting the unknown state p onto a coherent state 
a, one can perform position displacement, and rotation, R{9), operations on the 

unknown state and measure the overlap of the resulting state onto the ground state |0), 
which is a straightforward process in our experiment. Applying the displacement operator 

1 /2 

amounts to physically displacing the state a distance of x = 2x0 l^l, where xq = {2^) 
is the ground state width of a particle of mass m in a harmonic oscillator of frequency 
u. In a harmonic oscillator R {9) = exp [— 2a^a6'] can be implemented by letting the state 
evolve under Ti = huj [a) a + |) for a time t = 9/uj. How these operations are performed 
experimentally and how the ground state population of a state is measured will be described 
in section III. 

The measurement of the Wigner distribution is simplified by recognizing that for any 
symmetric non-degenerate potential and, in particular, in a harmonic oscillator, W (0, 0) = 
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n I-^ (^1 P 1^^) ~ ^ Yl'^=o (^1 P I"-)- values away from the origin the displace- 

ment and rotation operators can once again be utilized: 

^ oo 

W{x,p) = W{2xo^e{a),2po'^m{a)) = - V (n |a) , (2) 

n=0 

where po = and p{n\a) = {n\ D'' {\a\) R (6) pK^ (6) D {\a\) \n) is the probability 

of finding the particle to be in state n after the application of the rotation and displacement 
operators to the unknown state p. Therefore, the determination of the Wigner distribution 
is reduced to performing population measurements after applying displacement and rotation 
operators. This is similar to the reconstruction of the Husimi distribution except that the 
final population measurement requires the measurement of all states, (|0) , |1) , |2) , ...), not 
just the population of the ground state. 

III. IMPLEMENTATION 

A. Optical Lattice Specifications 

A 1-D optical lattice is formed by the counter-propagating components of two laser 
beams resulting in an intensity interference pattern of the form I (x) = Iq cos^ (fcxsin (|)), 
where k = ^ is the wave vector of the laser and 7 is the angle of intersection of the two 
beams. This standing light wave induces a light shift on the atoms resulting in a potential 
U (x) = I (x) where A is the detuning of the laser light from atomic resonance and 
Fand Ig are the natural line-width and saturation intensity of the atom respectively. Of 
importance to this experiment is the fact that the individual wells of an optical lattice can 
be approximated as harmonic oscillators. Therefore the theory of Husimi and Wigner distri- 
bution measurements in harmonic oscillators can be applied here. The oscillation frequency 
in each well 00, and lattice depth Uq, are related by the equation u = ^^sj^- We find uo 
directly by measuring the period of Ramsey fringes created by inducing oscillations with a 
spatial shift of the sinusoidal potential {2^. 

Our experiment starts with a cloud oi^^Rb atoms in a magneto optical trap (MOT) which 
are then cooled in an optical molasses to a temperature on the order of lOpK. The optical 
lattice is turned on during the MOT stage as cooling in the presence of the lattice increases 
the lattice loading efficiency. The two laser beams intersect at an angle of 7 = 49.6°, 
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resulting in a lattice vector of fc/, = ^sin(^) = f = 3.38 ■ 10^m~^, where A = 780nm 
is the wavelength of the lattice light and a = 0.93/im is the spatial period of the lattice. 
The lattice has a detuning of A ^ 27r ■ 25GHz from the F = 3 ^ F' = 4 D2 trapping 
line of ^^Rb so as to make the scattering rate negligible and to allow the sinusoidal lattice 
potential to be treated as conservative. The lattice is tailored to support 2 — 4 bound states, 
depending on the experiment, which means a depth of f/o — 10 — AOE^ where Ej. = 
is the scaled recoil energy in the lattice direction. The lattice is formed in the vertical 
direction so that atoms in unbound states will, in a time of < 10ms, leave the interaction 
region due to the pull of gravity. Each of the lattice beams passes through an acousto-optic 
modulator (AOM). By controlling the frequency and phase of the signal with which we drive 
each AOM independently we are able to control the position, velocity and acceleration of 
the lattice. Shifting the relative phase of the lattice beams by displaces the lattice (or, in 
the rest frame of the lattice, displaces the atoms) by a distance d = a0/27r. Spatial shifts, 
used for the displacement operators D (|a|), can be applied with a resolution on the order of 
Inm. Displacements as large as twice the lattice spatial period occur in a time of < 0.5/is. 
With typical oscillation frequencies of a;/27r ~ lO'^Hz, this time interval for a shift can be 
considered to be instantaneous. 



B. State Preparation 

The preparation of the states whose quasi-probability distributions we shall measure 
requires two steps. The first is to filter out the ground state. By lowering the intensity of 
the lattice beams over a time of 10ms until only one bound state is supported, and keeping 
it there for 5ms, all atoms in higher states become unbound and fall out of the interaction 
region due to gravity. Afterwards we raise the intensity back to its original level but are 
left with the ground state in each well, with a typical contamination of the first excited 
state of about 5 — 15%. A more detailed explanation of our ground state preparation and 
measurement process can be found in reference |3| and a similar technique is described in 
reference j^l- The implementation of the second part of the state preparation depends on 
what state we wish to prepare. 

In this paper we prepare and measure three states: a ground state, a near coherent state 
and a state with a population inversion. For the ground state no further action is necessary 
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as it was prepared in the filtering stage. A near coherent state, \j3) , is prepared by shifting 
the potential by 5a; = | = 0.155/im (or a 60° phase shift). The magnitude of [3 is related 

1 /2 

to the displacement by = 5x (^) = 0.88. For both the ground and coherent state a 
lattice depth of ST-E,. was chosen, supporting 4 bound states with an oscillation frequency of 
uj = 4:8.33kHz. The deviation of the near coherent state which we prepare from an actual 
coherent state is due to the finite depth of the lattice. With only 4 bound states only the first 
4 terms of the coherent state are present in the lattice. This does not significantly detract 
from the validity of the approximation, as the amplitude of the n^^ term in a harmonic 
oscillator coherent state is exp [- \pf /2] ^ |/?|", which when \p\ = 0.88 is quite small for 
n > 4. This state was allowed to rotate in the potential for a time t = 20/is, giving it a 
rotation of 6* = 0.97 radians. 

Creation of an inverted state begins with the preparation of the ground state, this time in 
a lattice containing only two bound states (u = 32.2kHz). Next, the potential is given a 60° 
phase shift (0.155/im), held there for 80/is and then shifted back to its original position. We 
have found that this process excites a large number of atoms into higher states^]. After 
waiting several milliseconds in order to let unbound atoms leave the interaction region we are 
left with what will be shown to be an incoherent mixture of ground and first-excited-state 
atoms in a ratio of roughly 3 to 7. 

C. State Measurement 

Here we describe the process used to determine the Husimi distribution for the ground 
state and coherent state. As per equation 1 the measurement of Q {\a\ ,9) takes three steps. 
First, we allow the atoms to undergo free evolution for a time t to let the state oscillate in 
the harmonic-like potential for a rotation of 9 = ut in phase space. We use a total of 27 
angles (or wait times, t) separated by 0.24 radians (or 5/is), spanning a range of 9e [0,27r] 
(5te [0, 130] /is). Second, a displacement of a; = 2xo |a| is applied to the lattice. This is 
the displacement operator, i5(|a|). A total of 19 different displacements are used, each 
separated by 25.8nm or a 10° phase shift of the potential for a total range of xe [0,465] nm 
or [0, 180] degrees. Along with the phase component, the total number of measurements is 
then 513 not including repeated measurements for statistical analysis and averaging. Next, 
a projection of this new state onto the ground state is performed. This first requires the 
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same process that was used to filter the ground state during the state preparation: the 
lattice laser intensity is lowered until only the ground state remains. After waiting ~ 20ms 
in order to let the unbound atoms become spatially separated from those still bound the 
atoms are illuminated with resonant light, and their fluorescence collected in a CCD camera. 
With this image the relative population of ground state atoms, as a fraction of total atom 
number, is measured. 

Figure 1 is a phase space diagram of the Husimi distribution for the coherent state. As 
expected, it is a Gaussian displaced from the origin by |/?| = 0.88 at an angle of roughly 
55°. Since a coherent state is merely a ground state displaced from the origin in phase 
space both distributions should have the same width. Figure 2 shows a cross section of 
the Husimi distribution for the ground state along with one of the coherent state. Fits 
of these Gaussians give rms widths of 144nm and 147nm and peaks of 0.267 and 0.257 
for the ground and coherent states respectively. Husimi distributions, by deflnition, are 
normalized. By integrating the curve shown in flgure 1 we flnd from j Q {\a \ ,6)\a\d\a\d6 = 



Q \2x' — — drdO = 1 that Xrms = 9Q.3nm. Here we find a surprising result. The 



we find Xrms = 96.3nm = 1.09-Xo. We believe that this discrepancy is due to inhomogeneities 
in the intensity of our lattice beams. We have separately measured the dephasing time T2 
to be on the order of 2 or 3 oscillation periods in our lattice implying Am ^ OAuj. 

As Xq = h/2muj, we expect the rms width of our atom clouds to be Xrms = {h/2m){\/uj) = 
Xq[1 + (Au/uj)^ + {Auj/ujy + . . .)]. For Au/u ~ 0.4, this means Xrms ~ 1.08x0, consistent 
with our observations. 

In addition, the peak of the Husimi distribution of a true ground state is K Our measured 
distribution has a peak value of 0.267 = 0.84/7r, obtained from a direct measurement of the 
ground-state population to be 84%. A harmonic-oscillator state with 16% population in the 
first excited state, when the ground-state width is taken to be Xrms as calculated above, has 
a Husimi distribution with a width of 141 nm, consistent with our measured value. 

Following equation 2 the Wigner distribution for the inverted state is measured by first 
applying the rotation and displacement operators R{0) and D (|a|), as in the Husimi recon- 
struction. For the rotation operators we sample 41 different rotation angles separated by 0.16 
radians spanning a range of 9e [0, 27r] equivalent to waiting times in the range 6te [0, 200] fis 
with a resolution of 5/is. We again use 19 different spatial shifts for the displacement op- 




ground state width of a harmonic oscillator of this frequency is Xq = 




88.0nm but 
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erators with a range of xe [0, 465] nm giving a total number of 779 measurements for the 
Wigner measurement. 

After the rotation and displacement operations are implemented, the state populations 
must be measured. Since the lattice holds two bound states any atoms excited to higher 
states by the rotation and displacement operators become unbound. After waiting 5ms in 
order to let these atoms become spatially separated we once again lower the lattice laser 
intensity until only one bound state remains. After waiting another 20ms in order to let these 
first-excited-state atoms leave the interaction region we then capture a fluorescence image. 
From this image we determine the ground and first-excited-state populations as well as how 
many atoms were in states n>2. Therefore we are able to construct the Wigner distribution 
with the following caveat: measurements of W (x, p) become increasingly uncertain for values 
for which the displacement operator D(|a;|) excites more atoms into n>2, since only the 
first two terms of equation 2 are known exactly. 

Figure 3 shows a cross section of our best estimate of the Wigner distribution: W' {x, 0) = 
^ (p(0 |q;) — p(1 \oc))- At the origin we see that the distribution is negative, a nonclassical 
signature which is characteristic of a population inversion. The inset shows the absolute 
upper and lower bounds, based on our data, of the Wigner distribution due to our lack of 
knowledge about the populations of states higher than n = 1. The upper bound is obtained 
by assuming that all atoms lost during the displacement and rotation operators were in 
the 2nd excited state; ^ (p (0 \a) — p{l\a) +p{n > 1 \a)), while the lower is constructed by 
assuming they were in the 3rd; ^ (p (0 \a) — p {l\a) — p (n > 1 \a)). Of course neither one 
of these extremes can be the case as they are not physical. It is known that the Wigner 
distribution for any state must go to zero as x,p oo. In addition, one can get a sense of 
the magnitude of the disagreement between W and W by investigating the normalisation 
of the three curves; J W {x,p) dxdp — 0.42, which is clearly not consistent with 1. The 
same normalisation integral, however, gives values of 10.32 and —9.49 for the upper- and 
lower-bound curves, respectively. Clearly, the real Wigner function lies between W and 
the upper bound, but much closer to the former. Theoretically, one expects a significant 
additional contribution from the second excited state for |a| ~ 1, but for larger values of 
|a| this should be cancelled out by higher populations in the third excited state. Near the 
origin, there was no population in states with n >2, so there is no uncertainty. 

Figure 4 shows a 3-D representation of the Wigner distribution measurement. Note the 
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rotational symmetry, which indicates a lack of coherence between energy eigenstates. Usually 
processes such as the displacement operator create coherent superpositions of states such as 
in the case with our presented Husimi distribution. The lack of any coherence properties 
here is due to the combination of dephasing processes inherent in our optical lattice 3| 
and the long delay between the population inversion creation and the time of measurement. 
Therefore the result of our state creation, as shown by figures 3 and 4 is an incoherent 
mixture of ground and 1st excited state atoms with most atoms being in the latter, as 
evidenced by the negative component of the Wigner distribution. 



IV. CONCLUSION 



Making use of a newly developed technique to measure state populations of the vibrational 
states of atoms in l-/im lattice wells, and our ability to perform arbitrary translations 
in phase space, we have reconstructed Husimi and Wigner phase space distributions for 
atoms in ground, coherent, and inverted states of oscillation in an optical lattice. This 
is the first complete quantum characterisation of the state of motion of atoms in such a 
system. The ground and coherent states are largely consistent with expectations, possessing 
essentially the same shape and width, although they indicate both some admixture of the 
excited state, which is understandable in light of our preparation procedure, and an apparent 
underestimate of the width of the ground-state wave function relative to the experimental 
measurement, which is probably due to inhomogeneities in the lattice beams. We observe a 
nonclassical signature in the Wigner function, reaching a maximum negative value of —0.12 
at the origin, in the case of the inverted distribution. 
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Figure 1: Phase space representation of the Husimi distribution of a coherent state 1 ,6) 
for \l3\ = 0.88, 9 = 0.97. From the contour plot in the inset one can see the displacement of 
the Gaussian from the origin. Axis units are in x/2xrms where Xrms — 96.3nm. 
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Figure 2: Cross sections of the Husimi distributions of the ground state (closed circles) 
and coherent state (open circles). As expected they have virtually the same width but 
the coherent state is displaced from the origin. In addition, the height and width, when 
compared to the true ground state of a harmonic oscillator of this frequency, show that 
there is some contamination of the state by first excited state atoms. 
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Figure 3: Cross section of the measured Wigner distribution. Negative values at the 
origin are indicative of a population inversion. Inset shows absolute upper (+) and lower 
(x) bounds of the distribution as explained in the text. Axis units are in x/2xrms where 
Xrms — llQAnm 




Figure 4: 3-D image of the measured Wigner distribution. Of note is the cylindrical 
symmetry showing a lack of coherence between the component energy eigenstates. 
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